A QUANTITATIVE INVESTIGATION OF THE POMERON 



E. Gotsman ('^)'[|E.M. Levin and U. Maor^'^)' [] 

^"•^ School of Physics and Astronomy 
Raymond and Beverly Sackler Faculty of Exact Sciences 
Tel Aviv University, Tel Aviv 

^''^ DESY, Hamburg, Germany 
and 

St. Petersbourg Nuclear Physics Institute, Gatchina, Russia 



Abstract 



A comparative investigation of various Pomeron models is carried out through 
the study of their predicted values of atot, B, and in high energy pp and pp 
scattering. Our results strongly support a picture of the Pomeron in which we 
have both moderate blackening and expansion of the p(p) - p amphtude in impact 
parameter space as a function of energy in the ISR-SSC domain. In particular, 
we obtain an excellent reproduction of the data with a hybrid eikonal model 
which combines the hard Lipatov-like QCD Pomeron with the old fashioned soft 
Pomeron and Regge terms. Our analysis shows that the additive quark model, 
at least in the naive form, is not compatible with the data. 
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At present there are two incongrous descriptions of the Pomeron. One is that 
based on the conventional phenomenology of soft hadronic scattering and the relevant 
calculations are carried out within the framework of Reggeon field theory . The 
soft Pomeron is constructed from multiperipheral hadronic (Regge) exchanges and has 
ap{0) =1. This is not compatible with the high energy experimental observation of 
rising hadronic cross sections. The old fashioned soft Pomeron was, thus replaced § 
by a soft supercritical Pomeron with an intercept ap{0) -1 > . This model provides 
a reasonable description of pp and pp scattering in and above the ISR energy range 
. It's principle assumption, that the supercritical Pomeron is a simple isolated pole 
in the complex angular plane, lacks field theoretic justification, and contradicts our 
experience with the multiperipheral model which fails to produce a pole with ap > 1 

The alternate description of the Pomeron is within the framework of QCD. Lipatov 
[^] has shown that perturbative QCD, in which a hard Pomeron is built out of multipe- 
ripheral high transverse momentum gluon exchanges, predicts a different Pomeron i.e. 
a series of poles in the complex j plane above unity at 1 < j < 1 + e . The resulting 
Lipatov Pomeron has a more complicated form jj^^ with a > 1 . A consequence of a 

pole with intercept a{0) is that atot oc s"'^^'^'^. Hence, both the soft supercritical and 
the hard QCD Pomeron predict a powerlike rise of the total hadronic cross section. 
This may continue until the unitarity bound gives rise to screening effects which satu- 
rates the growth of the total cross section |^ making atot ^ In^s, which is compatible 
with the the Froissart limit 0. 

In the following we shall examine a wide class of models of the Pomeron. Our goal 
is to compare the predictions of these models with the high energy experimental data 
presently available over the energy range 23 < ^/s < 1800 GeV. The purpose of 
such a comparative study is to attempt to discriminate between models with different 
high energy predictions based on data which are apparently below asymptotia. 

A quantitative investigation of the various Pomeron models can be readily carried 
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out in the impact parameter space. Our amplitudes are normalized such that 
7r|/(s,t)p ; atot = 47r Im f(s,0) where 

f{s,t) = ^Jdhe'^"a{b,s) (1) 

and 

a(.,6) = ^/dqe-^'^-V(s,t) (2) 
Zn J 

Hence we have atot = 2 / rfb Im a(s,b) and aei = J dh\a{s,b)\'^ . The corresponding 
forward slope and curvature parameters are 

^ Jdhb^a{s,b) 

2jdha{s,b)' ^ ' 

^_ 1 /dbb^a(g,b) 1 ^ dhb'^a{s,b) ^ 

32 jdha{s,b) 16 ' j dha{s,b) ' ^ ^ 



2 



where 

^d<y^ ^{Bt+Ct^+ ) 



Unitarity requires Im a(s,b)< 1. In order to satisfy this constraint it is convenient 
to express a(s,b) in terms of the complex eikonal function h) with 

a(s,6) = i[l - e'^(^'^)] (5) 

where Im x ^ 0- This ensures that unitarity is restored on summing all the eikonal 
multi-particle exchange amplitudes. At high energy, forward elastic scattering is es- 
sentially diffractive and therefore Rex is very small. We assume that Rex ~ 0, conse- 
quently the amplitude a(s,b) is purely imaginary and determined by the opaqueness 
Vt{s, b) = Imx ■ This assumption limits our study to the forward direction, where the 
bulk of the experimental data is present. 

We would like to note that considering only the pure imaginary Q{s,b) violates 
the crossing symmetry of the scattering amplitude. We will describe later how we 
restore this sjTumetry, assuming that the real part of the amplitude is small enough. 
One of our major objectives is to critically examine how much of the observed change 
in the forward high energy pp and pp scattering amplitudes can be associated with 
blackening as opposed to the expansion of Q{s, b) and a(s,b) with s. We note that a soft 
Pomeron trajectory is essentially fiat i.e. a' is either zero or very small. Accordingly, 
the growth of the cross section associated with the soft Pomeron is attributed mostly 
to a blackening of the b-space scattering amplitude. This property is common to 
many other models |^ |T^. The b-space characteristics of the Lipatov Pomeron are 
less clear. In the LLA one expects |12[ a moderate expansion with energy. We will 
check whether this result is compatible with data. 

Even though we choose an eikonal approach we do not attach any profound theo- 
retical importance to this choice. It is merely a method to restore s-channel unitarity, 
and it enables us to assess the natural scale of the shadowing corrections. These can 
manifest themselves either on the hadronic level with a typical scale Rh ~ 1/m or on 
the quark scale with a typical value of Rq ~ 0.1/m. We expect the data analysis to 
provide us with some information on this. We would like to mention here that the 
eikonal formula (5) allows us to incorparate in our fitting procedure new ideas, such 



as were discussed in refs. [0 , |jT8[ about the shadowing correction using different 
assumption about Q{s,b). 

To this end we compare several Pomeron models with three experimental quantities, 
which are particularly sensitive to the above, i.e. atot , B and . We note that more 
than (Tei itself, the ratio is sensitive to the hadron opacity. Over the ISR-Tevatron 
energy range this ratio was found ^ to be in the interval 0.15- 0.25, which is well 
below the saturation limit of 0.5 corresponding to a fully absorbing black disc |^. 

The above experimental quantities relate directly to the structure of the inelastic 
(multiparticle) high energy amplitude. Indeed, atot is driven by the inelastic amphtude, 
B is the mean radius of the partonic distribution in the impact parameter space and 
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characterizes the hadronic opaqueness which is very sensitive to the scale of the 
shadowing corrections. 

We examine the entire statistics of 74 data points |^ of atot, B and obtained 
from to pp and pp scattering in the ISR-Tevatron energy range. In our handhng of the 
data, we have used the pubhshed results as such, without any attempt to discriminate 
or average between different ISR data sets at the same energy, or re-analyse the raw 
data as recently advocated In particular, our initial supposition, that a(s,b) is 
imaginary, prevents us from calculating p = without additional assumptions. 

The data we have examined are characterized, in general, by relatively small errors. 
Consequently, one must carefully evaluate the quality of the fit to this data. In our 
study we have found that seemingly good reproductions of this data set, yield an 

2 

unexceptable ^ > 10. An additional fitting problem is associated with the fact that 
most of our data comes from the ISR, with the UA4 and Tevatron providing only 

2 

6 high energy points. As a result we had to check, that fits with an acceptable 
also had a reasonable high energy behaviour. We would also like to mention that 
the detailed high t structure of the elastic cross section is beyond the scope of our 
investigation. This corresponds to the very fine details of a{h, s) or ^2(6, s) at small b, 
a region which we do not pretend to investigate. 

We have attempted a number of generic forms for the Pomeron, with different 
degrees of complexity. As we are discussing both pp and pp scattering, we have also 
included in our parametrization , an odd Regge term to account for the pp and pp 
difference. We therefore consider an amplitude of the form 

a{s,h)=Ap{s,h)±AR{s,h) (6) 

Initially, we attempt an orthodox non-unitarized parametrization of the Pomeron. 
This may possibly be adequate, since the measured data is well below the unitarity 
limit, and there may be no need to use unitarized expressions in our analysis. At this 
level we have for the soft Pomeron and Reggeon 

F,(s,t) = Qe^o.*s"'W-iszn[^a,(t)] (7) 

where i = P or R. With linear trajectories aiit) = aj(0) + a'^ t, we can readily take 
the Fourier transform of Eq.(7) and obtain 

Ai{s,b) = 

c^^g-p^^^'t-^^^'] ■ {Rlstn[-aM + Z,] - -a[cos[-am + Z,]) (8) 



where 



= Rl + a[lns (9) 
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Accordingly, Ai{s,b) has four adjustable parameters : Ci, Rl^, ai{0) and a[. 

The first model that we have tested has a supercritical soft Pomeron, and we 
denote it as aj (s,b). Such a pomeron can be viewed as the consequence of the soft 
multihadron (multiparton) production at high energies. The best fit parameters and 

2 

the ^ for this model are summarized in the upper section of Table I. In our original 
fit we found that a'p ~ 0.3, ~ 1. and a_R(0) ~ 0.5. We have, therefore, fixed these 

2 2 

parameters at their traditional values to somewhat improve our ^ . The ^ = 2.04 
is misleading. The model, as such, provides a good description of the ISR data at the 
expense of a very poor reproduction of the UA4 and Tevatron data. This problem was 
traced to a low fitted value of e = ctp - 1. We have also tested a model in which e is 
fixed at 0.085 and denote it as an (s,b,e = 0.085). The quality of the fit, presented 
in the upper section of Table I, is considerably poorer, but yields a fair continuation 
from the low to the high energy domains. 

The third model we have tested is that based on the Lipatov-like hard Pomeron, 
where we set 

with 

= a4 + a5(/ns)°« (11) 

The Qi denote six adjustable parameters. The parametrization considered is 

aiii{s,b) = ALi{s,b)±AR{s,b) (12) 

This parametrization assumes the soft contribution to be entirely associated with the 
Regge sector whereas the hard contribution is given by the Lipatov-like Pomeron driven 
by hard multihadron (multiparton) production. 

The form we have choosen for Ali is very similar in it's energy dependence to the 
form suggested by Bourrely, Soffer and Wu 0. The two models are, however, funda- 
mentally different. In the Bourrely, Soffer, Wu model, A(s,b) has a fixed b distribution 
independent of s. Accordingly, the rise in atot is due only to the blackening of A(s,b). 
The Lipatov 1 amplitude that we have examined, e.g. Eq. (10), has a b distribution 
which is energy dependent, and hence R^^ is s dependent, see Eq. (11). 

2 

Our best fit with the corresponding ^ is summarized in the upper section of Table 
I. As can be seen, the hard Lipatov-like Pomeron provides an excellent fit, which is 
considerably better than the soft Pomeron descriptions. However, it is important to 
note that neither of these models are suitable parametrizations at very high energies 
as the power dependence of s will eventually violate unitarity. This is particularly true 
for the Lipatov Ali{s, b) where we have A^i^s, 6 = 0) > 1 for -^i > 62 TeV. 
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To overcome this problem we have also examined the above three Pomeron models 
in an eikonal parametrization, where we set 



nj{s,b)=Ap{s,b)±An{s,b) (13) 
nii{s, b, e = 0.085) = Ap{s, b, e = 0.085) ± Ar{s, b) (14) 
^lni{s,b) = ALi{s,b)±AR{s,b) (15) 

Our motivation for utilizing the eikonal expansion rests, not only on the need to uni- 
tarize our amplitudes and study the shadowing effects , but also on the knowledge 
1^ that the eikonal description provides a natural explanation for the changing log- 



arithmic slope of ^ with t. Some prudence must be exercised at this stage. Our 
investigation aims at understanding the general characteristics of the Pomeron. As 
is well known, neither of Eqs.(13 - 15) can reproduce the dip structure at ISR and 
its slow change to a shoulder as observed by UA4. However, these are exeedingly 
small phenomena where ^ at the second peak is 10~^ — 10~^ smaller than ^ in the 
forward direction. Accordingly, when viewed in 6 - space it corresponds to a very 
subtle substructure of VL{b, s) at small 6, which is completely outside the scope of our 
investigation. 
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Our best fits and the corresponding are presented in the middle section of Table 
I. We again note, that the Vlj fit agrees well with the ISR data , but fails to reproduce 
the high energy data, due to the the low value of e in this fit. The quality of VLui fit 
is considerbly worse than the one obtained without eikonalization. We have therefore 
attempted a hybrid model 

fi,y(s, h) = Aii(5, b) + Ap{s, b) ± Ar{s, b) (16) 

Namely, we associate the increase of atot with the hard QCD Lipatov-like Pomeron 
which is superimposed on a non- leading soft Pomeron with ap(0) = 1, plus a Regge 
contribution. As can be seen from the middle section of Table I, we obtain an excellent 
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fit with ^ = 0.96. This is demonstrated in Figs. 1-3 where we show the hybrid Vtiy 
fit compared to the relevant experimental data. The energy scale was choosen so as to 
include our predictions for LHC and SCO. The quality of our fit is even more impressive 
if we consider the fact that some of the ISR data points, at the same energy differ 
by more than one standard deviation. The one data point which is not in agreement 
with the fit is the value as measured by UA4. We call attention to the fact that 
an alternative analysis of the UA4 raw data [|] yields a ratio which is compatible with 
our fit. 

A basic conclusion that can be drawn from our analysis, with or without eikonal- 
ization, is that Giot growth is associated with both the expansion and the blackening 
of a(s,b) and f2(s, 6) with s. In the hybrid model we have 

^sic^=23'fe'=o)'^ ~ ^"^^ ^w\^/^2z) =1-48. An identical result was obtained by Chou 



and Yang [11 1 who found, in the framework of their model, the need to change both 



the normalization and the radius of their f2(s, 6) as a function of s. 
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The results of the Lipatov 1 sector of the hybrid model can be readily compared 
with some theoretical expectations 0, 12 1 in which 



with 



= ao^=exp ^ (17) 
ylns 



In the above one should distinguish between parameters that are determined as a 
direct consequence of perturbative QCD, and those which depend on the boundary 
conditions and the LLA used by the authors of references [Q, |T^. On the qualitative 
level our fit supports the Lipatov energy dependence ^^^"^^3 with both 02 and 03 positive 

as expected. We also obtain R\i{s) oc ^/Ins, in agreement with references |12|. On 
the other hand, we see that the actual numerical values of our fitted 02,^3,04 and 



05 differ from the numbers calculated by references |T^. This difference can be 
traced to the fact that the above calculations were made in the LLA with specific 
boundary conditions. Both may require some modifications. We find it encouraging 
that Rl^{y/s = 23) ~ 17 GeV~^, which corresponds to a radius of Ifm. It is interesting 
to examine both the low and the high energy limits of the hybrid model. We note that 
Lipatov 1 in the low energy limit predicts B values which are considerably smaller than 
the experimental data. This is compensated, however, by the non leading soft sector, 
producing an overall reproduction of the data. This is not suprising, as the hybrid 
model suggests a smooth, rather than an abrupt (threshold), continuation from soft to 

3 

hard physics. Below the black limit, Lipatov 1 produces a (Ins) 2 dependence of atot on 
s. This is less than the {InsY allowed by the Froissart limit |]5[. In order to investigate 
this point further, we have also examined an alternative Lipatov 2 parametrization in 
which of Eq. (11) is replaced by 

R?i2 = fl4 + a^Vlns + a^lns (19) 

and we have 



nv{s, b) = Al2{s, b) + Ap{s, b) ± Ar{s, b) (20) 

This model has an {InsY dependence built in, provided Og > 0. Our best fit with this 

2 

form is ^ = 1.02, is also given in the middle sector of Table I. We find Og = 0.10 
wheras 05 = 0.789, rather close to the LLA estimates. 

Both fljv and fly models produce a curvature C which changes sign at -^i ^ 1200 
GeV and becomes more negative with increasing s. This result which is compatible 
with other models, is just another indication that the data examined is far from asymp- 
totia. We note that the non-eikonalized models predict very small C values over the 
entire range 23 < ^/s < 1800 GeV. Unfortunately, a direct experimental measurement 
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of C at t = is not possible, and the extrapolation of C measured at some finite t to 
t=0, is not reliable. 

As a check of the parametrization obtained in fit Q/y ( where the parameters were 
determined from fitting data in the energy range 23 < < 1800 GeV ), we have also 
evaluated the predictions of the model for atot, ^ , and B for lower energies i.e. 5 
< < 23 GeV. The results are in very good agreement with the experimental data 
are shown in Figs. 1-4. 

As a further check on our assumption of taking Re f(s,t) ^ for small values of 
t, we evaluate the differential cross-section obtained from the contribution of Im f(s,t) 
only. Prom the equality Im f(s,0) = ^ , we take 

~ <L (21) 

^c/r*"" 167r ^ ^ 

so that 

|-j^<^L-exp(K + a^) (22) 

for evaluating the differential cross-section at small t values. In Fig. 5 we show how 
our best parametrization describes the small t - behaviour of the differential elastic 
cross section. We interprete this figure as support of the approximations used, and to 
show how close the physical picture we extract from this fitting, is to reality. 

We now discuss the real part of the scattering amplitude which is intimately con- 
nected with the crossing symmetry. Our basic assumption was that the real part is 
small. We now check this assumption and calculate the real part for our parametriza- 
tions. To do this, we separate the different contributions of positive and negative 
signature in the whole amplitude / (s, t) ( see Eq. (1)). In our case, this is easy to 
do, since the negative signature contribution is due only to reggeon exchange. Finally 

f{s,t) = Uis,t) ± f.{s,t) (23) 

for pp and pp respectively. 

1. The real part of the negative signature contribution is given by the usual formula, 
namely 

Ref^is. t) .TT , , , 

Imf.is,t) = ^^-(2"^— • (24) 

2. For the positive signature we calculated Ref+{s,t) as follows 

f{^(;'^ = Q2 = - coti^c^nt = 0)) . (25) 
/m/+(s,t = 0) ^2 + ^ ^ ^ 

where a^^{t — 0) was calculated by fitting /+(s,t = 0) as a power s°+^(*=°) in the 
wide range of energy. Fig. 6 shows the result of the above calculation for our the best 
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fit and illustrates that p = f^jf^f^^-^ is small enough to support the assumption that 
one can neglect the real part of the amplitude for small values of momentum transfer 
squared. One can also see that we get a good description of the experimental data 
that have not been included in the fitting procedure, except of the UA4 data point at 
a/s = 546 GeV, which has been a suject of much debate f^. ^ 

An alternative method for examining the various suggested forms for the Pomeron 



is based on the additive quark model |Tj] . In such an approach we introduce ^ two 
scales for the hadron: it's size with a radius of Rh ~ ^fm, and the size of a constituent 
quark with Rq ^ 0.1 — 0.2/m. Based on this slightly naive picture, we can then follow 
Chou and Yang |Ty] and write the scattering amplitude in the form 

f{s,t) = QG\t)fQQ{s,t) (26) 

where G'^{t) = (j^y with z/^ = 0.71 GeV^, and Jqq denotes the quark-quark scatter- 
ing amplitude. We now apply our assumptions for the Pomeron and Regge structures 
to the QQ amplitude. 

We assume a simplified static model [0, ^ where the QQ interaction is observed by 
static spectator quarks. In such a model, we define a'^'^{s, b) and Q'^'^{s, b) in analogy 
to Eq.(5), obtaining cr^t^,cr^'^, Bqq and Cqq. The hadronic properties are derived 
from: 

(Ttot = 9(t2,? (27) 
B = Bqq + ^ (28) 

C = Cqq + -^ (29) 



1/4 



From the relation 



^ / (ite^*+^*' (30) 
167r i ^ ' 



we can calculate the ratio We have examined this model in detail, utilizing the 
various parametrizations shown before, and were unable to fit the data. We consis- 

2 

tently obtained -jj- > 10, which is indicative of our inability to reconcile, atot , the 



slope B, and the curvature parameter C in such a model. It has been suggested 0, |T5 
that the QQ amplitudes approach the black disc limit. To test this, we have also 
examined a parametrization in which 

aQQ{s, 6) = (1 - ^) ■ [eia^lns + 04 - b]] + 



^ The knowledge of Imf{s,t) and p enables us to calculate the real part of the amplitude at any 
t utihzing ^ 

Ref{s,t) = p^^{tlmj{s,t)) . 

The above expression completes our restoration of the crossing symmetry in the case when p is small 
enough. 
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+ e[b - a^lns - aj ■ exp[-^^-^^^^^^^] ± ag^(s, h) (31) 

as 

2 

Here again, we obtained a > 10. Our consistent inability to obtain a reasonable 
fit , is a strong indication that the additive quark model, at least in the naive form we 
have studied, is not amenable to our problem. 

Finally, we have also examined the possibility that the increase of atot with in- 
creasing s is connected to the onset of minijets produced by parton-parton collisons 



at energies above the ISR region [jT6[. In such a model the b-space transform of G {t) 



is normalized by the hard parton-parton cross sections to which we also add a soft 



background. We disregard the fact that the various authors of reference ||T6[ actually 



calculated an inclusive rather than a total cross section, and examine the generic form 

It is obvious that in this type of model, the increase in atot with s is due to an increase 
of Vt through blackening rather than expansion, which will only start when a(s,b = 0) 
= 1. Our results for this parametrization are summarized at the bottom of Table I. 

2 

The = 3.7 we have obtained implies that such a model cannot accurately account 
for the observed data. We therefore disagree with the conclusions of Block, Halzen 
and Margolis who have normalized such a model to the Tevatron data in order to 
make predictions for energies in the LHC and SSC domain. Our analysis, which shows 
both a blackening and expansion of a(s,b) and fi(s,6), excludes the possibility that 
the onset of minijets is the exclusive reason for rise of atot with energy. 
We summarize with the following conclusions and remarks; 

1) Our parametrization denoted by fi/y provides a very good description of the ex- 
perimental data over the energy range 5 < ^ < 1800 GeV. 

2) Regardless of details, the data in the 23 < i/i < 1800 GeV range shows both mod- 
erate blackening and expansion of a(s,b). 

3) The data is compatible with a smooth transition from a soft to a hard Pomeron 
contribution which can account for the rise of atot with s. We find that the onset of 
hard minijet production cannot exclusively account for the rise of atot- Our analysis 
provides strong evidence that the Lipatov Pomeron is the proper way to describe the 
rise in atot which is due to hard processes, such as minijet production. 

4) The additive quark model, at least in it's naive form, fails completely to provide a 
simultaneous fit to the data. 

5) The various models we have discussed produce somewhat different values of atot 
, B and at \/s =40 TeV. These are summarized in Table II. We note that our 
predictions for atot are higher than the 121 mb predicted by reference |3] and the 125 
mb predicted by reference It is interesting to check our asymptotic predictions at 

^ After the completion of this work we learnt of a more sophisticated approach to the quark model 
(see refs.|Q that takes into account the interaction between more than one quark from each 
hadron. 
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the Planck scale s ~ 10'^'^ GeV^. We observe that, asymptotically, the various models 
compared in Table II have reached the black limit where = 0.5, and have converged 
to the same limiting cross section atot = 1010 mb. If we compare this value with the 
Froissart bound 

o-toi = ^In^— (33) 
/i^ So 

we find that /i ^ 3 GeV for < 23 GeV. 
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Table 1 : Values of parameter for different parameterizations of the Pomeron. 
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Figure captions 

Figure 1 atot in mb vs ^/s in GeV. The measured cross sections are from . The 
dashed (sohd) hne are the predictions for pp { pp ) for fit Qjv 

Figure 2 B, the measured and predicted nuclear slope parameter in units of (GeV)"'^ 
vs ^/s in GeV. Data from 0. Dashed (solid) curve as in Fig. 1. 

Figure 3 The ratio ^ vs. ^/s in GeV. Data from . Dashed (solid) curve 
as in Fig. 1. 

Figure 4 The predicted value of the curvature parameter C in units of [GeV)'"^ vs. 
^/s in GeV. Dashed (solid) curve as in Fig. 1. 

Figure 5 Elastic differential pp cross section at small -t , at representative en- 
ergies. Data from . 

Figure 6 Ratio of real to imaginary parts of forward elastic amplitude vs. i/i 
in GeV. Data from 0. Dashed (solid) curve as in Fig. 1. 
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